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1. 


INTRODUCTION  AND  SUMMARY 


This  paper  is  concerned  with  the  problem  of  best  estimation,  i.e., 
minimum  variance  unbiased  estimation  (mvue),  of  certain  functions  of  the 
parameter  0 and  the  truncation  points  of  a power  series  distribution 
(PSD)  truncated  on  the  left  at  several  known  or  unknown  truncation 
points.  More  specifically,  given  k independent  random  samples  X , 

ij 

J = 1)  i = 1)  •••5  k,  the  ith  sample  from  the  PSD  with  prob- 

ability function 


a. (x)  0X 

(1.1)  p. (x)  = P[X. . = x]  = — 

1 1J  fi(8,ri) 


r.+l, 

l 


N. 

l 


where 


N. 

l 

(1*2)  f^COjr^)  = £ a (x)  0X,  (N.  finite  or  infinite)  , 

t=r.  1 

l 


it  is  desired  to  estimate  some  function  of  0,  including  p. (x)  itself; 

mvue  are  also  provided  for  the  truncation  points  r , . . . , r . 

i k 

Estimability  conditions  for  certain  functions  of  the  parameter  0 
of  a one-parameter  PSD  with  probability  function 

(1.3)  P[X  = x]  = , x e T (0  > 0) 

where  TCIq,  I the  set  of  non-negative  integers,  were  given  by  Patil 
[l963]»  Thus  it  was  shown  that,  on  the  basis  of  a random  sample  of  size 
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n >_  1 from  (l„3),  0 is  estimable,  i.e. s it  has  an  unbiased  estimator 

if,  and  only  if, 


(1.10 


n[T]  = Smin(n[T])  5 


where 


n[T]  = T +T2+0„„+T^.  with  T4  = T,  (i  = 1,  . » , , n) 

and  S denotes  the  tail  set  of  a,  that  is, 
a 

S = {x  : x e I and  x > a}  „ 
a o — 

In  particular,  (l.U)  implies  that  0 is  not  estimable  whenever  T is 
a finite  set,  e.g. , 0 = p/q.  in  the  binomial  case.  On  the  other  hand, 

0 is  estimable  if 


T - S , r e I 

where  r denotes  the  truncation  point,  assumed  known.  This  situation 

includes  the  Poisson,  the  logarithmic  series  and  the  negative  binomial 

distributions  truncated  on  the  left  at  a known  point.  In  fact,  Roy  and 

ot 

Mitra  [1957]  provide  the  mvue  of  9 for  any  integer  a > 0 in  the 
form 

•g  _ A(x-a,n,r ) 
a A(x,n,r) 


(1.5) 


where 


(1.6) 


A(x,n,r)  = ^ a(x  ) • ••  a(x^) 


the  summation  extending  over  all  n-tuples  (x^ , . „ „ , x ) of  integers 

x.  > r with  xn  + . . . + x = x. 
i — 1 n 

The  corresponding  case  when  the  truncation  point  r is  unknown  was 

studied  by  Josh!  [1972],  who  showed  that  on  the  basis  of  a random  sample 

x,  , . ...  x from  (1.3)  the  mvue  of  r is 
I n 


(1.7) 


/w  \ _ A(x,n,y+l) 

r y,X  y A(x,n,y)  - A(x,n,y+l) 


where 

x = x + .0.  + xr  and  y = min(x^,  . x^)  . 

Interestingly  enough,  though  9 is  not  estimable  in  the  case  of  a 
finite-range  distribution,  it  was  shown  that  the  probability  function 
itself  (hence  also  the  corresponding  distribution  function)  of  a PSD 
truncated  on  the  left  is  estimable,  both  when  the  truncation  point  r 
is  known,  Patil  [1963],  or  unknown,  Joshi  [1972]. 

The  main  difficulty,  reflected  even  in  the  treatment  of  some  special 
cases  of  truncation  away  from  zero  (r=l),  Tate  and  Goen  [1958], 
Cacoullos  [1961],  Patil  and  Bildikar  [1966],  Cacoullos  and  Charalambides 
[1972a],  lies  in  deriving  an  explicit  or  computationally  attractive 
expression  for  the  distribution  of  the  sufficient  statistics,  the  sum  of 
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the  observations;  equivalently,  to  get  around  the  evaluation  of  expres- 
sions such  as  A(x,n,r)  in  (l„6)„ 

The  above  difficulty  was  resolved  by  employing  exponential  gener- 
ating function  (egf)  techniques,  which  emerge  quite  naturally  in  the 
treatment  of  the  aforementioned  distribution  problem,  the  n-fold  convolu- 
tion of  a PSD  truncated  on  the  left  at  r(r  _>  l)„  Indeed,  Char alamb ides 
[1972],  [1973a] , [l973bl  and  Cacoullos  and  Charalambides  [1972b]  show  how 
the  egf  technique  provides  a unified  approach  for  a solution  to  the 
problem  of  mvue  of  certain  functions  of  0 in  (1.3),  both  when  r is 
known  or  unknowns  At  the  same  time,  the  treatment  of  truncated  versions 
of  classical  families  of  PSD!s  leads  to  natural  extensions  of  certain 
kinds  of  numbers;  the  Poisson  family  is  associated  with  generalized 
Stirling  numbers  of  the  second  kind,  whereas  the  logarithmic  series 
distribution  with  Stirling  numbers  of  the  first  kind  (the  well  known 
Stirling  numbers  correspond  to  the  simple  case  of  truncation  away  from 
zero,  i.e,  , r-~l);  finally,  the  binomial  and  negative  binomial  cases  led 

to  the  introduction,  Cacoullos  and  Charalambides  [1972a],  and  extension, 
Charalambides  [1972],  [1973a],  of  a new  kind  of  numbers,  called  C-numbers, 
The  more  general  situation  of  truncation  at  several  known  truncation 
points,  that  is,  (1,1)  with  a^s.x)  — a(x),  i = 1,  t , k,  motivated  the 
introduction  of  certain  multiparameter  Stirling  and  C-numbers,  Cacoullos 
[1973a]  and  [1973b],  The  mvue  of  0a  and  the  probability  function  itself 
were  given  in  terms  of  these  multiparameter  numbers , Cacoullos  [1973c], 

Here  we  consider  the  relevant  distribution  (Section  2)  and  estimation 
problems  (Section  4)  in  the  more  general  case  of  (l,l);  especially,  when 
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the  truncation  points  r , . ..,  r are  unknown,  a. (x)  = a(x)  and 
N.  = “,  i = 1,  . ..,  k.  When  p = (r  , . . . , r ) is  known,  the  mvue  of 
p.(x)  is  also  considered.  In  Section  3 we  review  the  multiparameter 
Stirling  numbers  and,  in  particular,  the  C-numbers  emerging  out  of  the 
more  general  set  up  of  different  a^(x)  for  the  binomial  and  negative 
binomial  distributions.  Some  of  the  results  of  this  Section  are  based 
on  properties  of  the  multiparameter  Stirling  and  C-numbers,  studied  by 
Cacoullos  [1973b],  In  Section  5*  we  exhibit  several  properties  of  condi- 
tional distributions  for  PSD's  and  indicate  how  they  can  be  employed  to 
provide  confidence  intervals  for  certain  functions,  in  particular  products 
and  quotients  of  the  parameters.  Applications,  including  confidence 
interval  estimation  of  attribute- failure  reliability  models  for  multi- 
component  systems,  are  exhibited  in  Section  6. 


2.  CONVOLUTIONS  OF  LEFT-TRUNCATED  PSD'S 


Let  X.  „ . i'  = l,  ...,  k j=l, 

X J 

with  p . f o ' s 


n^  be  independently  distributed 


(2.1)  p.(x)  = P[X 


x]  = 


ai  (x) 

f.^r.) 


ex, 


x = r„,  r . +1 , ....  N . 

i i l 

j = 1 , . . . , ni 


where 


is  independent  of  the  parameter  0 and  the  series  function 


N. 

f‘i(05ri)  = l ai(x)6X 
x=r . 

l 
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converges  in  some  interval  0 < 0 < ; N may  be  finite  or  infinite 

(i  = 1,  o„o9  k)„  When  the  truncation  points  r_^  are  known,  the  k-sample 
sum 


k 


n. 

OL 


X 


- I I 

i=l  «j=l 


ij 


is  a complete  sufficient  statistic  for  0,  since  a PSD  belongs  to  the 
exponential  family  of  distributions . It  is  easily  shown  that 
(cf.  Cacoullos  and  Charalambides , [1972 ]) 


Proposition  2.1.  The  p.f.  of  X is  given  by 


[2.2)  P[X=x]  = 


c(x;v,p)  Qx 
Sv(03p)  x! 


x = m,  m+1 , . . , , M , 


where  we  set  V = (n^,  , » , , n^ ) , p = (r-^,  „ , , a r^). 


(2„3)  m = n^ r^  + „ * » + 


\rk  = M = n„Ni  + + nkNk  » 


and  g (0,p)  is  the  egf  of  the  numbers  c(x;V,p)  given  by 


(2A) 


M x k 

gv(tso)  = l c(x;p,v)  — = n ~ [f.tt^r,)] 
x=m  " ~ ” i=l  i ' 


n. 

l 


moreover,  the  numbers  c(x;V,p)  have  the  representation 


(2,5) 


c (x;V,p ) = 


xl 


n ^ 

r 


k 

T I n 


n„ 

i 


II  a. (x. „ ) 
n.  I ° . , „ , l ij 

k x 1=1  j=l 


6 


the  summation  extending  over  all  n-tuples  (n  = n + . . . + n ) of 

^ -L  K 

integers  x. . such  that 


x„ „ > r „ 

ij  — i 


I 


i=l 


n. 

1 


0=1 


ij 


x 


Suppose  now  that,  in  addition  to  0,  the  truncation  vector  p is 
unknown.  Then  by  the  results  of  Fraser  [1952],  a complete  sufficient 
statistic  for  (0,p)  is  (X,Z)  = (X,Z  , „ „ , , Z ) where 

~ ~ _L 


Z . = min  {X.  . } i = 1 , 


!<j<ni 


For  our  purposes , we  require  the  following 


Definition  2.1.  Let  G(u)  be  a real- valued  function  with  domain  the 
set  of  lattice  points  u = (u^,  u^)  with  integer  coordinates. 

The  kth-order  partial  difference  of  G(u),  to  be  denoted  by  A^,  is 
defined  by 


A G(uj  = (A  Au  o , o A ) G(u) 
u ~ u.  k-1  u. 

k 1 


where  A denotes  the  usual  forward  difference  operating  on  u. , 
u.  i 

l 

that  is , 


(2.6) 


A^  G(u)  = G(u+e3 ) - G(u) 
i 


denoting  the  ith  coordinate  axis  unit  vector  (0 , „ . <,  90,sL,Q , „ . . ,0) 
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The  following  lemma , of  interest  in  itself,  will  be  used  in  deriving 
the  p,f,  of  the  sufficient  statistic  (X,Z). 

Lemma  2.1.  Let  Q(u)  denote  the  tail  probability  function  of  an  integer- 
valued random  vector  U = (U^,  . U^),  that  is, 

Q(u)  = P[U  > ux,  . . . , Uk  > uk]  = P[U  > u]  . 

Then 

(2.7)  P[U=u]  = (-1  )kA  Q(u)  . 

~ ~ u ~ 

Proof:  For  k = 2 , we  have 

P[U1=u1,U2!Su2]  = Q(u1,u2)  - QC^+ljUg)  - Q(ulSu2+l)  + Q(  ^+1,^+1)  , 

which,  in  view  of  (2,6),  gives  (2,7)  - For  k > 2 the  proof  is  analogous 
(cf.  the  kth  order  partial  difference  A^F(u)  of  a k-variate  continuous 
distribution  function  F(u)  which  is  used  in  the  definition  of  the  density 
at  u by  letting  h = (h  , h ) -*  0). 

•v  ~ K 

¥e  can  now  state 

Theorem  2,1,  The  distribution  of  the  sufficient  statistic  (X,Z)  is  a PSD 

(_l)kAzc(x;z,v)  „x 

(2.8)  P[X=x  ,Z=z  ] = ~ (Q^p)"  (v,y)  <x<M,  l IP 

where  g^(0,p)  and  c(x;z,p)  are  defined  by  (2.3)  and  (2A). 
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Proof:  Setting 


(2.9) 


Q(z;x)  = P[Z>_z,X=x]  , 


we  obtain  by  Lemma  2.1 


(2.10) 


P[X=x,Z=z ] = (-1)KAzQ(z;x) 


However , we  have 


Q(z;x)  = P [ Z>_ z , X=x]  = P[X.  ,>z.  , X=x,  i=l,.  . . ,k,  j=l,.  . . ,n. ] 
~ ~ ~ ij — l i 


, n. 
k i 


n [f  (0,r  )] 
i=l 


— I n n a. (x. . ) 


i Z i=l  j=l 


where  the  summation  extends  over  all  n-tuples  of  integers  x^  with 


. n. 

k i 


x_„  . _>  z^  and  \ \ x.  . = x 


ij 


i=l  j=l 


ij 


Thus , by  (2.3)  and  (2.4),  Q(z;x)  can  be  written  as 


(2.11) 


Q(z ;x) 


c(x,z,v)  qX 
SV(0,P)  x! 


Hence,  as  A operates  only  on  z,  we  have 


A Q(z;x)  = 


A c(x;z,v) 

5 ~ ~ ex 


DV 


l,p)  x!  * 


and  (2.8)  follows  from  (2.10), 
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3. 


MULTIPARAMETER  STIRLING  AND  C-NUMBERS 


It  was  already  mentioned  in  Section  1 that  the  present  approach  to 
the  distribution  theory  relating  to  mvue  problems  of  a left-truncated  PSD 
was  motivated  by  the  results  of  the  simple  case  of  a Poisson  distribution 
truncated  away  from  zero,  which  leads  in  a natural  way  to  the  well-known 
Stirling  numbers  of  the  second  kind  S(x,n),  with  egf 


g(t)  = £ S(x,n) 

x=n 


In  the  present  notation 9 this  corresponds  to  the  Poisson  series  function 

f (X,l)  5 f(X,l)  = eX-l  (k=l)  . 


The  general  case  of  truncation  of  a PSD  distribution  on  the  left  at 
several  points  r^,  , r^  motivated  (Caeoullos,  1973 ) the  definition  of 

multiparameter  numbers;  their  egf  is  always  of  the  general  form  (2.4). 

Of  special  interest  are  the  following  cases: 

a)  Poisson:  Here  0 is  the  Poisson  parameter  A and  the  PSD  series 


functions  are 

f,  (A 

i 


r.-l 

) = f(X,r  ) = eX  - l 

j=0 


jl 


i = 1. 


k . 


The  corresponding  egf  defines  the  multiparameter  Stirling  numbers  of  the 

second  kind  S (x;p,v);  for  k = 1 they  reduce  to  the  generalized 
K. 

Stirling  numbers  of  the  second  kind  S(x,r^,n^)  considered  by  Charalambides 
(1972);  for  k = 1 and  = 2 they  give  the  associated  Stirling  numbers 
of  the  second  kind  (see  Riordan  [1958],  p.  77). 
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Recurrence  relations  can  be  obtained  by  using  the  difference-differen- 
tial equation  satisfied  by  the  egf  G (t,p),  say,  of  the  numbers  S(x;p,v) 


du  V 


G.jC'tsp)  = nGv(t,p)  + l 


r . -1 
l 


A <v1)! 


(t,p) 


where,  as  before,  n = nn  + „ „ . + n^.  In  the  sequel,  we  will  require  the 
(x,v)-wise  recurrence: 

k / x 

(3.1)  S(x+l;p,v)  = nS(x;p,v)  + £ I 

i=l  \ i- 

b ) Logarithmic  series:  Here  we  have 


S(x-ri+l;p,V-ei) 


r .-1 
J A J 

f(0,r  ) = - log(l-e)  - l i = 1,  k (r.  >1) 

J=1  J 1 


where  0 < 0 < 1.  The  corresponding  egf,  h (t.,p),  say,  defines  the  multi- 
parameter  signless  Stirling  numbers  of  the  first  kind,  |s(x;p,v)|;  the 
numbers  s(x;p,v),  the  multiparameter  analogues  of  the  usual  Stirling 
numbers  of  the  first  kind  s(x;n),  are  generated  by  the  basic  series  func- 
tion log(l+0) , instead  of  - log(l-0),  and  their  egf  is 


hv(t,p) 


n„ 

°°  tX  1 

I s(x;p,v)  it  = n 

x~m  ~ ” i=l 


1 

n.  ! 

1 


log(l+t ) 


r.-l 

1 

I 

j=l 


(-l)J-1 


The  corresponding  expansions  (representations)  of  S(x;p,v)  and 
|e(x;p,v)j  (cfc  (2„i))  show  that 

s(x;p,v)  = (-l)n_X|s;(x;p,v)  | . 
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For  k = 1,  = 1,  s(x,l,n  ) and  js(x,l,n^)j  reduce  to  the  usual 

Stirling  and  signless  (positive)  Stirling  numbers  of  the  first  kind, 
respectively „ 

Recurrence  relations  for  s(x;p,v)  and  )s(x;p,v)|  can  be  obtained 

from  the  corresponding  difference-differential  equations  satisfied  by 
* 

their  egf's  h^(t,p)  and  h_ (t ,p) , respectively: 

, * k r,-l  r.-l  # 

(1+t)  ~ hv(t,p)  = l (-1)  1 t 1 hv_e  (t,p) 

i=l  ~ ~i 

k r . -1 

d-t)  ^-hv(t,p)  = l t 1 hv_e  (t,p) 

i=l  ~ ~i 


Thus  we  obtain  the  recurrences : 


k r4  -1 

a(x+l;p,v)  + x 3 (x;p9v)  = l (-1)  "t  (x)  _1(x-ri+l  ;p  9 V-e± ) , 

i=l  i 


(3  = 2) 


|s(x+l;p9v)|  -x|s(x,p,v)|  + j (x)  |s(x-r.+l;v-e. ,p) 

i=i  i 


c ) Binomial : Now  we  have 


f „ (0 ,r . ) = (1+0 , 


r.-l  Is.' 
s „ 1 / i 

) 1 “ I 


= 1,  . . o , k 0 = p/q. 


J=0  \ J 


(the  s,  denoting  positive  integers).  The  corresponding  egf  H (t,p) 
i V ~ 

defines  the  multiparameter  C-numbers  C(x;p,V,a)  with  a = (s  ) 

/v>  ~ ~ X.  K 


The  special  simple  case  of  k = 1 with  s^  = s , n^  = n,  r^  = 1 , 


motivated  the  definition  of  the  C-numbers  CX  (Cacoullos  and 

s ,n 
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Char alamb ides  [1972a]);  in  the  present  notation,  we  have 


(3.2a) 


C(x;l ,n,s ) = sn  C*  n = ^ (-1  T"'3  (J)  ( sj  ) 


This  explicit  expression  was  obtained  by  elementary  occupancy-model 
arguments  in  obtaining  the  n-fold  convolution  of  a binomial  distribution 
(with  parameters  s and  p = 0/(l+0))  truncated  away  from  zero  (also 
referred  to  as  positive  binomial). 

Recurrence  relations  for  the  C-numbers  can  be  obtained  from  the 
difference-differential  equation  (of  Cacoullos,  1973 ) 


(3.3)  (1+t) 


_d_ 

dt 


H (t,p,a)  = (v,a ) H (t ,p,a] 

v ~ ~ ~ v ~ ~ 


i=l 


(S.  )r 

l l 

(r.-l)  i 


r .-1 

t 1 H (t ,p,a) 

V— c:  „ ~ ~ 

- -1 


Thus  an  (x,v)-wise  recurrence  is 


(3.4) 


C(x+l;p,v,a)  = [(a,v)-x]  C(x;p,v5a) 

k / 

’ \ 


i=l 


x 


1 / 1 


It  should  be  noted  that 


C(x;p,V,a)  =0  if  x < m — (p,v) 


and  the  C-numbers  are  positive  integers  for  m < x < (v,a). 
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d)  Negative  Binomial s Here  we  have 


__  g 2.  I — 3 \ 

f. (e,r. ) = (1-6)  1 - l I 1 (-1)3  9J  1 = 1,  ....  k 

1 1 3=0  \ 3 I 

(&  > 0 and  0 = p;  for  the  Pascal  distribution  the  are  integers). 

The  corresponding  egf  defines  the  signless  C-numbers 

|c(x;p,V,-a)  = (-1)  C(x;p,V,-a);  the  numbers  C(x;p,V,-a)  are  defined 

by  the  egf  H^(u;p,-cf),  that  is,  as  in  e)  with  a replaced  by  -0, 
so  that  x has  an  infinite  range  of  values 


x >_  (p,v)  . 

It  can  be  shown,  e.g.  by  using  the  analogues  of  (3.3)  and  (3 A),  that 
(3.5)  |c(x+l;p,Vs-<7)|  = [(c,v)  + x]  !c(x;p,v-e.  ,-o)  j 


k / x \ 

+ y 


r . -1 

\ 1 / 


(s.+r.-l)  C(x-r , +1 :p ,V-e » ,-d) 
i i r. 1 l ~ ~ ~ 

i 


It  should  be  observed  that  in  the  special  ease  k = 1 , r^  = 1,  n^  = n, 
S1  = S’ 


C(x;l»n,-s) 


vx  n 


n' 


' ^ l (_l)n“J  [ J (sj+x-1)^  , 


ns 


x 


since  in  (3.2a)  (sj)  must  be  replaced  by  (-sj ) = (-l)X(sj+x-l)  . 

X X X 

For  additional  recurrence  relations  and  combinatorial  interpretations 
of  the  multiparameter  Stirling  and  C-numbers,  we  refer  to  Cacoullos  (1973b) 


- ih  - 


k.  MINIMUM  VARIANCE  UNBIASED  ESTIMATION 

For  the  construction  of  the  mvue  of  a parametric  function  4>(0,p) 

A c . 

it  suffices  to  find  an  unbiased  estimator  cp(x,z ) which  is  a function  of 
the  sufficient  statistic  (X,Z),  whose  distribution  we  considered  in 
Section  2.  Extending  the  results  of  Patil  [1963]  and  Joshi  [1972] 
concerning  the  case  k = 1,  or,  by  a close  examination  of  the  condition 
of  unbiasedness,  as  obtained  by  using  the  p.f.  (2.8), 

I I (-!)k  4>(x,z)  A c(x;z,v)  ~ = g (0,p)  <f) ( 6 , p ) 
x>m  z>p  t ~ ~ x!  Z 

we  get  the  following 


Theorem  U.l.  On  the  basis  of  the  k independent  samples  X.  „ , i=l,...,k 

j=l,.,.,ni  from  (2.1),  a parametric  function  (0,p)  has  a mvue  if,  and 
only  if,  for  every  p = (r1,...,r  ) and  every  0,  (a)  the  function 

<f>(09p)  gv(0,p)  admits  a power  series  expansion  in  0,  say 


Xk.2) 


4>(e ,p)  g (e,p)  = l c (x;p,v) 


xeT 


X! 


and  (b)  for  each  p the  index  set  T of  cj)(0,p)  g^(0,p)  is  a subset  of 
the  index  set  of  the  series  function  g (0,p),  i.e.  (see  (2.2)  and  (2.3)), 

(^•3)  T d{m,  m+1,  . . . , M>  . 

A 

If  <j>(x,z)  exists,  then  it  is  given  by 


(b.b) 


<j)(x,z) 


* 

AZC  (x;z,v) 
A b(x;z,v) 

Zj  'V-  ~ 


- 15  - 


Proof:  We  have  from  (4.1),  (4c 2)  and  (2.4)  for  each  p 


l l 0(x,z)  AzC(x;z,v)  6X  = l C (x;p,v)  0X  , 

x>m  z> p ~ ~ xeT 


which  will  hold  for  all  0 if,  and  only  if, 


(4.5) 


l (-1)  <j)(x,z)  A 0(x;z  ,v)  = C (x;p,v) 

*>p  t 


Applying  the  same  argument  involved  in  (2.9)  and  (2.10)  to  the  function 


Q (p)  = I (-1)  4>(x,z)  A^C (x;z  ,v) 

z-2p 


we  get  from  (4.5) 


* * 

A Q (p)  = 4>(x,p)  A c(x;p,v)  = A C (x;p,v)  . 

~ jj  ~ ~ M ~ ~ 


Since  this  holds  for  every  p,  (4.4)  readily  follows. 
From  the  preceding  theorem,  we  have 


Corollary  4.1.  The  truncation  point  p = (r  ),  as  well  as  . 

~ _L  K. 

p = (r^,..c,rk),  a > 0 integer,  is  always  mvu  edfc&mahle ; the  mvue  of 

(a)  . 

p is 


0i 

^ A [z  c(x;z,v)] 

(4,6)  P (x,z)  = — — 

-a  ~ A..c(x;z,v) 


A/ . nC (x;z+e  ,v)  a-1  . , n 
za  + _ LXJ ~ .~x  ~ y fa] 

' \ \ L 


i Azc(x;z,v)  „-Q  !■ 


i = 1, . o o ,k 
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where  , denotes  the  (k-l)-order  partial  difference  operator  n Az 

Ati  J 

It  should  be  noted  that  the  second  expression  in  (4.4)  gives  the 

"correction"  term  for  the  corresponding  maximum  likelihood  estimate  za 

i 

Ql 

of  r.  to  yield  an  unbiased  estimate  of  r „ 
i i 

Corollary  4,2  The  probability  functions  p^r  +j)  in  (l.l)  are  mvu 
estimable  and  the  mvue  of  p.(r.+j  ) is  given  by 


(4.7)  Pi(ri+j ;x,z)  = ai(zi+j> 


Az[ (x)z  +jc(x-zi-j ,z j V— e^ ) ] 


n^Az  c(x;z,v) 


j>0 

5 

i=l , . . „ ,k 


hence  the  corresponding  distribution  functions  are  also  mvu  estimable. 

As  in  the  case  of  known  truncation  points,  it  is  observed  that  condi 
tion  (4.4)  is  not  satisfied  if  all  the  N in  (l.l)  are  finite  and  the 
parameter  to  be  estimated  is  0 However,  we  have  the  following 


Co roil ary  4.3.  For  any  integer  a > 0,  0a  is  mvu  estimable  if  at  least 
one  of  the  N in  (l.l)  is  infinite;  that  is,  at  least  one  of  the 
distributions  p (x)  has  range  equal  to  the  tail  set  of  r. . Then  the 

1 i 

QL 

mvue  of  0 is  given  by 


(4.8) 


0 (x,z)  = (x) 
a 'a 


Azc(x-a;z,v) 
A c(x;z,v) 


for  y _>  p,  x _>  m . 


In  particular,  the  mvue  of  0 is 


(4.9) 


0 ( x , z ) 


A c (x-a;z,v) 
z ~ ~ 

Ac (x;z,v) 

Z)  -V 
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a 

For  comparison  purposes,  we  give  the  corresponding  estimates  of  0 
and  p (x)  obtained  by  Cacoullos  [1973c]  when  the  truncation  points 
r , ...jr^  are  known.  They  are 


(4.10) 


c(x-a;p,v) 

b(x;p,v) 


x > m + a 


(4.11)  pi(ri+j;x)  = ai(ri+j! 


(x) 

d (x-r . - j ; p ,V-e . ) 1 r.+j 


e(x;p,v) 


n. 

l 


x > m + j 

J > o 


for  i = 1 , , . . ,k  . 

It  is  interesting  to  observe  that,  the  role  of  c(x;p,v)  when  p is 

known  is  played  by  (-1)  A c(x;z,v)  when  p is  estimated.  The  mvue  of 

Z ~ v 

the  estimable  parametric  functions  in  specific  situations  are  obtained 
from  the  preceding  general  expressions  by  replacing  the  numbers  c(x;z,v) 
by  their  analogues  in  each  case:  in  the  Poisson  case,  c(x;z,v)  = S(x;z,v) 
(the  multiparameter  Stirling  number  of  the  second  kind);  in  the  logarithmic 
series  distribution,  c(x,;z,v)  = |s(x;z,v)j  (the  multiparameter  signless 
Stirling  numbers  of  the  first  kind);  In  the  binomial  case,  c(x;z,v)  = 
C(x;z,v,a)  (the  multiparamet er  C-numbers);  in  the  negative  binomial  case, 
c(x;z,v)  = | G (x;z ,v ,-c; ) j (the  multiparameter  signless  C-numbers). 

In  certain  cases,  the  mvue  obtain  a more  illuminating  form  by  using 
some  of  the  recurrence  relations  given  in  Section  3.  Thus  the  mvue  of  the 
Poisson  parameter  X can  be  written,  by  using  the  recurrence  (3.1),  as 
follows : 
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A(x) 


x 

n 


k 

x-1  \ 

1 - y 

L 

r . -1  / 

i=l 

l / 

S(x-r. ;p,v-e. )' 
x ~ ~ ~ i 

S(x;p,v) 


= c X (x) 
1 o 


when  p is  known,  and 


A(x,z)  = 


x 


k 

lx-!  \ 

1-1 

r.-ll 

i=l 

l i / 

£(x;z,v) 


= C2XQ^  » 


when  p is  estimated,  where  we  set  g(x;z,v)  = (-l)kA  S(x;z,v)  and 

~ ~ z ~ ~ 

/s  ~ 

Xq(x)  = x/n  = x denotes  the  usual  mvue  of  X under  no  truncation,  i.e., 
p = (0,0,,.., 0);  c and  may  he  regarded  as  the  correction  factors 

hy  which  Aq(x)  must  he  multiplied  to  give  the  corresponding  estimate  of 
X under  truncation  on  the  left  at  a known  or  unknown  point  p, 
respectively.  It  should  he  noted  that  0<c2<l.  Similarly,  in 

the  negative  binomial  case  the  mvue  of  0 (the  usual  p)  can  he  written, 
hy  using  the  recurrence  (3=5) 3 as  follows: 


9(x)  = 


x 


(a,v)+x-l 


k 

1 - I 

tx-l  \ 

i=l 

l / 

| C(x-r  ;p,V-e  ,-a) 

(S.+r.-l)  - 1 1 ~ 

ii  r . 

i 


C(x;p,V,-a) 


when  p is  known,  and  a similar  expression  with  |c(x;p,v,-a) | replaced 
by  A |C(x;z,V,-a) | when  p is  estimated.  In  both  cases,  0 takes  the 
form 


x 

c (a,v)+x-l 


A 

where  0q  is  the  usual  mvue  of  0 in  the  absence  of  truncation  and  c , 
o <_  c <_  1,  the  corresponding  correction  factor. 
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5.  SOME  CONDITIONAL-DISTRIBUTION  PROPERTIES  OF  A PSD 

In  the  present  section,  we  state  certain  properties  of  conditional 
distributions  of  PSD’s;  some  of  these  properties  have  been  used  for 
confidence  interval  estimation  of  the  reliability  of  multicomponent 
attribute-failure  models  in  certain  special  cases  such  as  the  negative 
binomial  (Hwang  and  Buehler  [1973])  and  the  Poisson  distributions  (Harris 
[l9Tl]);  apparently,  however,  the  general  underlying  property  for  PSD's 
has  not  been  explicitly  stated!  The  following  property  is  basically  a 
consequence  of  the  exponential  structure  of  a PSD, 


Proposition  5,1,  Let  X1 , „ » , ,Xg  and  Y^ , „ . , ,Y^  be  independently 
distributed  each  according  to  a PSD  and  let 


P F X „ -X  ] = 

1 


a.  (x) 

1 ' ' x 

A_.  (X,  } Ai 


x e S, 


(5^1) 


1 1 

b.(y) 


p[Yj=yi  ' l^UTT  “3 

<j  <3 


y e T 


j 


i = 1. 


j = 1. 


s ,, 


, k , 


where  the  sets  S.  and  T.  are  subsets  of  the  set  I of  non-negative 
1 .1  o & 


integers  and  each  of  the  series  functions 


.(V-  1 


xeS, 

1 


(x)X^ 


B (y  ) = l b (y)y^ 

J j y£T^  J J 


converges  in  some  interval;  Q < A„  < L„  , 0 < y„  < M„  (i=l,.„„,s  , 

11  J J 

j=l,...,k). 
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Define 


(5.2) 


u.  = x.  - xn 
111 

X ““  2 g e o o j S ^ 

D - <U2>. 

V.  = Y,  + xn 
J J 1 

J “l5ooo?k  £ 

y » (v1„ 

Then  the  conditional  distribution  of  X 


given  U = (u  , „ . . ,u  ) = u and 
X ■’v'  ^ s 


V = (v  , . . » ,v  ) = v is  a one-parameter  PSD  with  parameter 
~ X K. 


s / k 

e = n a . / n y „ 
i=i  x/j-a' 


and  probability  function 


(5.3) 


p0(x;u,v) 


4(x;u,v) 
D(0 ;u,v) 


with  series  function 


D(0;u,v)  = J d(x;u,v) 
xeT 


d(x;u,v)  = a.  (x)  II  a„(u„+x)  II  b„(v„-x) 
i=2  Jl=1  j j 


where  the  range  T of  x is  determined  by  the  ranges  S, , T and  the 

x j 

values  of  the  conditioning  variables  u and  v, 

The  proof  is  straightforward  and  is  therefore  omitted.  It  is  observed 
that  the  conditional  distribution  (5.3 ) depends  only  on  the  parameter  0, 
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not  on  the  individual  A^!s  and  ’ s 0 The  disappearance  of  all  nuisance 
parameters  is  guaranteed  by  the  general  result  of  Lehmann  and  Scheffe  (see 
Lehmann,  1959,  Ch.  2,  Lemma  8)  concerning  exponential  families.  It  follows, 
furthermore,  that  the  theory  of  uniformly  most  powerful  unbiased  tests  and 
corresponding  uniformly  most  accurate  unbiased  confidence  intervals 
(Lehmann,  1959,  Ch.  h)  can  be  applied  for  inferences  about  0, 

Another  conditional  distribution  arises  in  relation  to  making 
inferences  about  0 in  the  presence  of  unknown  truncation  points  playing 
the  role  of  nuisance  parameters.  Since  in  applications,  especially  in 
reliability  of  attribute-failure  models  (Hwang  and  Buehler  [1973],  Harris 
[1971]), 


we  are  going  to  consider  conditional  distributions  given  only  U,  i=e. 

'V 

with  conditioning  variables  differences  of  X^,!s;  in  the  preceding 
Proposition  5„1,  this  would  lead  to  a PSD  with  p.f,  (of  given  U=u) 

d(x;u) 

(5.M  ■ Dtqry 

where 

S 

d(x;u)  = a,  (x)  II  a„(u.+x)  , D(A;u)  = £ d(x;u)A'  . 
x i=2  1 1 x 

i=l,..o,s,  be  independently 


Proposition  $.2.  Let  X„ „ 
distributed  with  p.f.'s 


a.  (x) 

P[X„  „=x]  = — -7^- r XX 

10  WV  1 


x = r, ,r.+l, „ . . ,N. 
i x 5 i 

i = 


Then  the  statistics 


X.  = l X 
i , 


i=l 


. . W.  = min{X„ ,} 

ij  i , ij 

J 


i '»  )S 


define  a complete  sufficient  statistic  for  the  parameter  set 

{0lS...,0k,  rlS...,rk}  and  the  conditional  distribution  of  X given 

— Y = s " ° « ivg ) and  U = u (U  defined  as  in  (5.2))  is  a one— parameter 

PSD  with  parameter  A = A ,,..,A  and  p.f. 

i-  s 


(5.5) 


p(x;u,w)  = - 


c (x;u,w)  ..  x 


A_ 

C ( A ;u,w) 


where  max[W  , max  (W  .-U )]  < x<  min{N  SN  -u  . . ,N.  -u  },  and 
± 2<l<s  11  ± d d . s s 


w 

c (x;u,v)  = A c (x;w,  ,n.  ) II  A c.  (x+u.  :w.  .n, 

W,  1 1 1 . „ W , 1 11] 


'1  - - - i=2  -i  A A x 1 (Vx)!  * 


iX 


C (A;u,w)  = l C (x;u,w)  yy  . 
x 


Proof.  We  have 


P[X  =x,U=u,W=w]  = P[X  =x,W  -w  ] n P[X. =u. +x,W.  =w. ] 

~ ~ ~ ~ x j.,  x 3_— 2 11  ll 


since  the  pairs  (X1 jW^ ),..., (X^ ,Wp ) are  independently  distributed.  Now 
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applying  Theorem  2.1  with  k=l  for  each  (X^,W^),  we  obtain 


P[X1=x,U=u,W=w]  = (-IT 


A c,  (x;wn  5n,  ) .x  A c-  u,+x;w.  ,n.  ) , i 

w-.  1 1 1 An  s w,  x i 5 i 5 i X . 

1 1 „ x 1 


S(x1,r1)  xli"2 


gn  (\+^)l 

i 


with  A operating  only  on  w.  (i=l,...,s).  Similarly,  by  Theorem  2.1, 
Wo  jL 

i 

the  joint  distribution  of  U and  W is 


k ui+k 

Ac  (k;w_  ,n.  ) 8..  s Ac.  (u.+k;w  n,  ) 0. 

p[u=usw=w]  . I {-1)S  1.  1 !--4  n -i-i, — 

k S1(Xl*I'1)  1=2  giUl’rl>  (Vk)l 


so  that  (5=5)  follows. 


6.  SOME  APPLICATIONS 

Here  we  indicate  some  applications  of  the  results  given  in  the 
preceding  sections.  The  first  kind  of  applications  concern  the  mvu  esti- 
mation of  parameters  involved  in  a PSD  truncated  on  the  left.  This  calls 
for  the  theory  developed  in  Sections  2,  3,  and  especially,  4.  The  second 
kind  relates  to  testing  hypotheses  and  constructing  confidence  intervals 
for  products  and/or  quotients  of  a number  of  parameters,  each  associated 
with  a distribution,  usually,  of  the  same  family.  This  was  already  briefly 
discussed  in  the  preceding  section. 

As  regards  the  first  kind  of  applications , there  are  many  practical 
situations  in  which  sampling  is  naturally  restricted  to  truncated  PSD’s, 
especially  the  classical  ones:  Poisson,  binomial  and  negative  binomial, 

with  typical  truncation  away  from  zero.  For  example,  a multiply  truncated 
Poisson  may  arise,  e.g.,  in  estimating  the  accident  rate  X on  the  basis 
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of  reports  from  several  sources  (localities,  factories,  etc.)  where  source 

i provides  relevant  information  only  if  at  least  certain  number  r of 

i 

accidents  occur;  r may  be  fixed  and  known  or  unspecified.  Other 
examples  in  which  left-truncated  Poisson,  as  well  as  binomial  or  negative 
binomial  distributions,  may  arise  can  be  found  in  the  statistical  litera- 
ture and  some  of  the  references  given  in  this  article. 

We  give  an  application  in  the  area  of  reliability  of  multi-component 
systems.  It  is  assumed  that  each  of  the  m,  say,  components  of  a complex 
system  (a  missile,  an  aircraft,  a computer,  etc.)  can  only  either  perform 
or  fail  so  that  the  system  reliability  in  this  so-called  attribute  failure 
model  can  be  expressed  as  a function  of  Bernoulli  parameters.  Consider 
the  reliability  of  a "k-out-of-m  system"  (see,  e,g.,  Birnbaum  et  al.  [1961]) 
which  consists  of  m components  and  will  function  properly  if  at  least  k 
out  of  m components  function, where  1 £ k £ m;  k = 1 corresponds  to  a 
parallel  system  and  k = m to  a series  system.  The  reliability  of  such  a 
system  with  m independent  identical  components  is  given  by  the  tail 
binomial  probability  function 


(6.1) 


m 


R 


l “ P^l-P) 

i=k\-"  / 


k 


m-r 


where  p is  the  reliability  of  each  component.  Suppose  observations  (of 
the  number  of  failing  components)  are  available  for  such  k-out-of-m 
systems;  then  the  mvu  estimation  of  R provides  an  application  of  the 
theory  of  Section  4 (see  (4,11)).  It  should  be  observed  that  if  we  were 


simply  interested  in  a confidence  interval  for  R = R (p),  then  it  would 


- 25  - 


be  possible  to  convert  a confidence  interval  for  p to  a corresponding 

one  for  R (p),  which  is  an  increasing  function  of  p. 

K. 

Moreover  9 in  certain  situations  it  may  be  assumed  that  data  for  such 

k-out-of-m  systems  are  available  only  under  the  condition  that  at  least 

r,  l_<r<k,  of  the  m components  function  in  all:. instances . Then  the 

corresponding  (conditional)  reliability  of  the  system  is  given  in  terms 

of  R by 
K 

* R, 

R = — 1 < r < k < m . 

r R — — 

r 

Thus  in  effect  we  arrive  at  a truncated  binomial  model,  where  we  sample 
from  a left-truncated  probability  function 

(m)  x,_  >m-x 

P[X=x]  = — — E-4--~£" . x «=>  . 

K 

r 

* 

Then  in  estimating  we  can  use  the  mvue  based  on  (lt.ll)  or  (i+.7), 

depending  on  whether  r is  known  or  unknown. 

With  respect  to  tests  of  hypotheses  and  confidence  intervals  for 

products  of  parameters  of  PSD's,  we  recall  the  application  of  Harris 

[1971]  to  products  of  Poisson  parameters;  in  an  attribute  failure  model 

it  is  assumed  that-  component  i fails  with  probability  p^ , i = l,,,,,k 

and  hence  the  reliability  of  a parallel  or  series  system  can  be  studied 

in  terms  of  the  product  p Pg.o.p  ; if  component  i is  tested  separately 

n^  times,  then  for  "large"  n_.  and  p_^  "small",  the  problem  can  be 

reformulated  in  terms  of  the  product  X = XnX_o,.X,  where  X.  = n.p.  is 

12k  1 11 
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the  parameter  of  the  ith  Poisson  distribution  which  approximates  the  ith 
binomial  (n^p^),  i = l,...,m.  The  theory  of  Section  5 for  constructing 
uniformly  most  powerful  unbiased  tests  and  corresponding  uniformly  most 
accurate  unbiased  confidence  intervals  applies  to  such  products  of  Poisson 
parameters  even  when  the  underlying  Poisson  distributions  are  left- 
truncated  at  different  known  or  unknown  truncation  points.  When  the 
truncation  points  r^,...,r  are  known  Proposition  5.1  applies,  and  when 
they  are  unknown  Proposition  5.2  applies;  note,  however,  that  the  latter 
case  requires  taking  more  than  one  observation  from  each  of  the  k Poisson 
distributions . 
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